Evolution of a Bose-condensed gas under variations of the confining potential 
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We discuss the dynamic properties of a trapped Bose-condensed gas under variations of the confining field and 
find analytical scaling solutions for the evolving coherent state (condensate). We further discuss the characteristic 
features and the depletion of this coherent state. 
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The recent successful experiments on Bosc-Einstein 
condensation (BEC) in trapped ultra-cold alkali atom 
gases Jl]-||] open a unique possibility to investigate dy- 
namic properties of a Bose-condensed phase. Of par- 
ticular interest is the response of the system to time- 
dependent variations of the confining field. In an in- 
teracting Bose-condensed gas these properties are non- 
trivial. For example, if initially almost all trapped atoms 
are in the condensate (T = 0), then under adiabatically 
slow change of the trapping potential they remain in the 
condensate which now corresponds to the ground state 
of the system in the instantaneous trapping field. At 
the other extreme, a fast change of the potential from 
the initial to final shape brings the system to an excited 
superpositional (ES) state, where the admixture of the 
(final) ground state can be small. Then, even assuming 
complete isolation of the system from the environment, 
there is the question of how the correlation properties 
change. Especially interesting are those responsible for 
the reduction of the probability of inelastic processes due 
to the presence of a Bose condensate |4j . 

Another question concerns trapped gases with nega- 
tive scattering length. The trapping field stabilizes the 
condensate provided the spacing between adjacent trap 
levels exceeds the interparticle interaction . Will this 
hold in the ES state or will the system "collapse" ? 

Of principal importance is the evolution of a com- 
pletely isolated many-body system which proves to be 
in the ES state. Does the system undergo stochastiza- 
tion and imitate relaxation behavior, at least for a large 
number of particles and large interaction between them 
compared to the level spacing in the potential well? This 
question is related to the well-known problem of the ap- 
pearence of irreversibility in a quantum system with a 
large number of particles. 

To answer most of the above questions we first con- 
sider the evolution of a Bose condensate in a parabolic 
trapping potential V(r) = mto 2 r 2 /2 (m is the atom mass) 
with frequency u(t) varing from loq to u>\ . We use ideas of 
the analysis of the quantum motion of a particle in a har- 
monic oscillator with time-dependent frequency (see Q|) 
and find the solution of the time-dependent non-linear 
Schrodinger equation for the evolving coherent state. In 



certain cases our analytical results can be compared with 
numerical calculations, first performed in ref. ||. We fur- 
ther analyze the characteristic features of the evolving 
coherent state and discuss the problem of relaxation and 
the loss of coherence. 

Let us consider a Bose gas with a fixed number of par- 
ticles N in a symmetric harmonic potential with time- 
dependent frequency uj(t). We assume a pair interaction 
potential between atoms of the form U(r) — U8(r). In a 
3-d gas U — Anh 2 a/m, where a is the scattering length 
and TO the atom mass. The Schrodinger equation for the 
Heiscnbcrg field operator of atoms, \l/(r,i), reads 
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The field operator can be represented as a sum of the 
above-condensate part and the condensate wavefunction 
*o =< $ >, which is a c-number: = + * (see 
||). The equation for \& (r, t), obtained by averaging 
both sides in Eq. ([!]), in the mean field approach has the 
form: 



* +L>[|%| 2 + 2n']* , (2) 



where n'(r, t) =< W>. The mean field equation of 
motion for ^'(r,t) follows from Eqs. (P and (^) : 
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In Eqs. (||) and (||), due to the condition 7r|a| 3 <C 1 in is the 
gas density), we omitted the terms containing anomalous 
averages < >. 

Frequency variations change the time and distance 
scales in Eqs. (||) and (||). Let us formally introduce 
new operators taking this into account (p = r/b(t)): 

*o(r,t) = [&(t)r d/2 Xo(p,T(tj)exp (i*(r,t)), 

#'(r, t) = [b(t)]- d ^ 2 x'(p, T{t)) exp (i*(r, f)), (4) 

where d is the dimension of the system and the phase 

<S>(r,t) = (mr 2 /2h){b(t)/b(t)}. (5) 

The equations of motion for the operators xo and x' m 
new coordinate (p) and time (r) variables take the form 
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where n'(p,r) =< X* (P> T )x'{p> r ) >■ Note that the 
choice of phase in the form (j^) leads to the cancellation 
of terms proportional to V P x' and V P xo- It is important 
that the phase is the same for both *5>o an d 

Let us first consider d = 2. We choose r(t) and b(t) 
such that they are governed by the equations: 



r(t) = J dt'/b 2 (t'), (8) 

6(t)+^(t)fe(t)- Wo 2 /6 3 W, (9) 

where loq = uj{— oo) is the initial frequency. Then Eqs. (|^) 
and (Q) are reduced to a universal form 
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Xo + ^[|xo| 2 + 2n'] X o, (10) 
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Eqs. ([h]) and ([tl]) are universal in the sense that in 
the variables p, r the problem is reduced to an inter- 
acting Bose gas in a harmonic well with constant fre- 
quency. Once we find the solution in the initial poten- 
tial well, we know the answer at any t. One should 
only solve the simple equation (g) with initial conditions 
b(— oo) = oo) = 0. For example, the expression for 

^o(r,t) reads 



r 
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where p is the initial chemical potential and 5 , o(r) is the 
stationary solution for the condensate wavefunction at 
t— > — oo. Using the Bogolyubov transformation, general- 
ized for an inhomogeneous case (see, e.g., [||), on the ba- 
sis of Eqs. (|l0|), ( [TT| ) we can describe the evolution of the 
spectrum and wavefunctions of elementary excitations. 

Neglecting the excitations, Eq.([l2]) describes the 2-d 
(radial) evolution of 'J'o m long samples (axial frequency 
is much smaller than the radial one) under variations of 
the radial frequency, as to a first approximation one may 
omit the dependence of "fo on the axial coordinate. 

If the non-linear interaction terms in Eqs. (g) and (|) 
can be omitted the above universal scaling takes place for 
any d. This requires at least a small ratio of interparticle 
interaction to the level spacing in the initial potential: 
rj=n U/?iui <^.l, where 710=1^0(0, -oo)| 2 (see below). 

In the 3-d case for arbitrary 77 there is no universality. 
On the other hand, in the limit 7/ 3> 1, neglecting the 
excitations, Eq.(g) can be again reduced to a universal 
form. In this case the kinetic energy term in Eq. (p) is 



comparatively small and can be omitted (see below). 
Then, introducing, instead of (||) and (g), the equations 
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r(t) = J dt'/b 3 (t'), 
b{t)+u 2 {t)b{t)=J 2 /b i {t), 
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For t— > —00 Eq.(|l5|) has a solution xo(^) =Xo( r ) ex P( — *M*)> 
where x is given by the well known expression |10|O]: 
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and zero otherwise. As follows from Eq. (|l5|) , for any t 
the condensate wavefunction has the form 
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where b(t) is governed by Eq.(jT^) with oo)=l, &(— oo)=0. 

Eqs. (jl2|) and dl7|) conserve the norm /|W (r, t)\ 2 d d r = 
No, where Nq is the initial number of particles in the 
condensate. Universal solutions of Eqs. ( |l0| ) and ( JTT| ) 
conserve the norm J (|\E'o( r i t)\ 2 + n'jr, t)) drr = N. 

We should emphasize that Eqs. (|l^) and (|l7|) describe 
a coherent evolution of ^q. Generally speaking, it is very 
different from the condensate wavefunction correspond- 
ing to the ground state of the system in the potential well 
with an instantaneous value of uj or the final value oji, 
even if ui(t) returns to the initial frequency ujq. 

The time dependence of "to and normal excitations for 
d = 2 (or at any d for 77 1, neglecting the non- linear in- 
teraction terms) is determined by the solution of Eq.(M. 
The latter can be found in a general case (see, e.g., M). 
Note that the classical equation for a harmonic oscillator 
with time-dependent frequency: 



£(t) + w 2 (t)£(t) = 0, 



(18) 



leads to Eq.(||), if one sets £(t) oc bit) exp(±itOaT(t)). 
With £(t) = exp(idJo^) for t — > —00 and t interpreted as 
a "spatial coordinate" , Eq. (|l8|) is equivalent to the one- 
dimensional Schrodinger equation for the reflection of a 
particle with "energy" Wq from the "potential" uJ 1 — w 2 (t). 
At times where to is already constant (u>i) we have 
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where R is the reflection coefficient and 6 the phase. 
For slowly changing frequency (on a time scale r ^> 
the coefficient R is exponentially small, and b = 



y/ojo/tui. In this case the initial condensate is adiabati- 
cally transformed to the ground state of the system in the 
final trapping field, without oscillations. If the condition 
^o.iTo 3> 1 is not valid, for at least one of the frequencies, 
the scaling parameter 6(f) will oscillate with a constant 
amplitude given by Eq.(19). There is no damping of the 
oscillations of the condensate density, unless relaxation 
is included (see below). 

The instantaneous size of the evolving condensate is 
related to the initial size by ro(f) = r b(t) (r = r (— oo)). 
In the case of abrupt change of the frequency we have 



R= (lu -wi)7(w + wi) 2 , 5 = 0, 



(21) 



and the function b(t) oscillates from 1 to loq/uji. For 
luq 3> wi there is a large expansion of the condensate and 
then compression to the initial shape. At times t <C uj^ 1 
the expansion is practically free, and Eqs. (|l9|), ( pp| ) yield 

b(t) = (1 + uj 2 t 2 ) 1/2 , r(t) = (1/wo) arctan(w £)- (22) 



If u>i — the compression does not occur. We have an 
expanding condensate described by Eq.(22) at any t. 

At times t ^> ujq 1 the characteristic velocity of free ex- 
pansion v = r (t) = idor . As follows from Eq.(16), 
for rj 3> 1 the initial size of the condensate ro = 
(2/i/mwp) 1 / 2 , and we have vq = •v/2/z/rn. Since in this 
case fi = uqU , the velocity is determined by the interac- 
tion between particles. For rj <C 1 the initial size ro ~ la, 
where Iq = (h/mujo) 1 ^ 2 is the amplitude of zero-point os- 
cillations in the initial potential, and Vq ~ ^jTiLu$/m is 
much smaller than for r\ 3> 1. 



The same picture holds for the 3-d evolution of a con- 
densate in the opposite limit, where initially rj 3> 1, al- 
though we should use the scaling transformation follow- 
ing from Eqs. ([l3|), (|lj). The latter has no analog in 
the quantum theory of scattering and should be solved 
directly. We again obtain a periodic function b(t), and 
oscillations of the condensate density will be determined 
by Eq.(|l7|) which was derived assuming a large ratio of 
interparticle interaction to the kinetic energy term (IK 
ratio) in Eq.(||). For r\ ^> 1 the IK ratio varies as rj 2 /b(t) 
and definitely remains large if the frequency increases 
(b(t) < 1). As follows from Eq.([l4|), fast decrease of the 
frequency to u>i -C <^o leads to large oscillations with 
frmax ~ \/2/3tL>o/wi. The solution ( |17| ) will be valid if 
w i ^ ^o/v 2 - ^ times t <C ^f 1 , where the expansion is 



velocity of expansion vq w y / 4/x/3m. For uj\ -C too/v 2 
the IK ratio can become small and the scaling (0), in- 
stead of (14), should be used. Then the profile of the 
condensate density will change, but the parameter b(t) 
determining the characteristic size of the condensate will 
be very close to that following from Eq.(|l4|). 

Let us now describe the evolution of a 3-d Bose- 
condensed gas with negative scattering length. The ini- 
tial condensate will be stabilized by the trapping field if 
tj«l, i.e., the IK ratio is small |3]||. The prime sta- 
bilization factors are the presence of the gap ~ hujo for 
one-particle excitations and the existence of a large en- 
ergy barrier for quantum fluctuations leading to collapse 
H . Both are related to small values of rj. In the case of 
radial evolution of long samples the IK ratio in Eq. (|Io| ) 
remains constant. Hence, the evolving condensate with 
a < and initially small r\ will be equally stable with 
respect to collapse as the initial condensate. For the 3-d 
evolution with r\ -C 1 the IK ratio in Eq.(||) varies as 
i]/b(t) and decrease of the frequency (expansion) makes 
the evolving condensate even more stable. On the other 
hand, with large fast increase of the frequency (compres- 
sion) to u>i > uo/r], or adiabatic increase to lo\ > ujq/t] 2 , 
the parameter rj(t)/b(t) strongly increases and becomes 
of order unity. This can lead to instability of the conden- 
sate with respect to collapse. The principal difference of 
the uniform 3-d compression from the radial compression 
of long samples is attractive for comparative experiments. 

The evolving coherent state described by the wave- 
function ^o(r,f) is an ES state. For sufficiently large 
and fast change of the frequency the admixture of the fi- 
nal ground state in ^o is very small, which raises two 
questions: What happens with correlations character- 
istic for the static condensate in the absence of irre- 
versible processes, and how fast is the depletion of the 
evolving coherent state. Analyzing the first question 
we consider correlations responsible for the reduction of 
the probability of inelastic processes due to the pres- 
ence of the condensate. The event rate of an m-body 
inelastic process in a homogeneous gas v m — a m Z m Q, 
where a m is the rate constant, Q the system volume, and 
Z m =< [*t(0,0)] m [*(0,0)] m > the local density corre- 
lator P). For three-body recombination we have m = 3, 
and for spin-dipole relaxation m = 2. In the absence of 
condensate Z m = m!(rT)™ 1 , where n is the average particle 
density. At T — > in a stationary condensate the density 
fluctuations are suppressed and Z m = (n) m . Hence, v m 
in the condensate decreases by a factor ml M. 

In the spatially inhomogeneous evolving Bose- 
condensed gas, generalizing the above expression for 
the event rate we have v m (t) = a m J d 3 rZ m (r,t). 
The structure of the field operators in the correlator 
Z m (r,t) is determined by the scaling transformation 
([|). If almost all atoms are in the coherent state 
*o(r,*), then both Eq.Q and Eq.@ lead to v m (t) 



free, Eq.Q yields b(t) w yflfiujQt (t 3> uj 1 ), and the v m (—oo) (Q(—oo) /£l(t)) m , where v(— oo) is the inelastic 
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rate in the initial static condensate and the appearence of 
the quantity (0(— oo) / 'f^i))" 1-1 is a trivial consequence of 
changing the system volume: fl(t) = £l(-oo)b d (t) . This 
result shows that coherent evolution retains the effect of 
partial suppression of inelastic processes, characteristic 
for the static condensate. Since any loss of coherence 
will lead to an increasing inelastic rate, there is an inter- 
esting possibility to study the depletion of the evolving 
condensate through the measurement of the rates of in- 
trinsic or light-induced inelastic collisional processes. 

In an isolated system the relaxation of the evolving 
coherent state will be accompanied by the appearence 
of an effective temperature. Assuming zero initial tem- 
perature and wi<woj a complete depletion of the con- 
densate would lead to the effective temperature T c f ~ /i 
for 77 ^> 1 and T c [ ~ foiuo for 77 <C 1. In the 3-d case 
the BEC transition temperature in the final potential, 
T c ps hu>iN^ 3 p2[ . For 77 <C 1 the condition of complete 
depletion of the condensate, T c /T c f ~ (lui/lu )N^ 3 <1, is 
rather strong although there is an upper bound for the 
number of particles: N^^lo/a. For rp$>l, with /i = r]hujQ 
and 77 ~ (Noa/lo) 2 / 5 H, we obtain a weaker condition 
which can easily be fulfilled: 

T c /T ci ~ (wi/Vo)^ 73 * («i/wo)(io/a) a/8 JV ~ 1/lfi £ 1. 

The question of relaxation and the loss of coherence 
in a quantum system with a large number of particles 
has several non-trivial aspects, especially with regard to 
dynamic evolution of a completely isolated system. In the 
latter case we can discuss the imitation of stochastization 
and the relaxation picture, although there are a number 
of reasons for the real relaxation. Thus far our analysis 
has assumed the mean field approach. The relaxation 
can only appear beyond this approach: One should at 
least consider the exact Hamiltonian and include in the 
analysis the interaction terms proportional to U 2 . 

Fast decrease of the frequency to uj\ <C 0J0 brings the 
system to a high ES state of the final potential. In 
the case of a large interaction (i; > 1) the character- 
istic energy of expansion is ~ fj,. For maximum expan- 
sion (b « Wo/wi) the evolution is almost quasistation- 
ary and the size ro(t) is close to the size of stationary 
states with energy ~ \x in the final potential. The corre- 
sponding one-particle density of states g(p) w [i 2 /2(huji) 3 
(d = 3) will be very large, ensuring that the spacing be- 
tween adjacent levels Se ~ /i(?iwi//j) 3 is much smaller 
than the interparticle interaction at maximum expansion 
nominll ~ fj,(oji/u}o) 3 . The many-particle density of states 
grows exponentially, with the exponent depending on the 
number of particles N. Under these conditions even a 
small external influence can lead to "mixing" of states 
and insert irreversibility. 

The characteristic relaxation time in this case will be 
determined by the collisional time r c ^ [ncrv(n)]~\ where 
a = 8na 2 is the elastic cross section, v((J.) the particle ve- 



locity at energy //, and n~nomin- Strictly speaking, r c 
represents the minimum relaxation time and it would be 
interesting to find possibilities to observe a larger time 
of relaxation. One may assume that even in the absence 
of any external influence r c will be responsible for the 
formation of the state which imitates relaxation. The 
imitation can be promoted by the deviation of the exter- 
nal field from harmonicity or spherical symmetry. 

The characteristic time of dynamic evolution is deter- 
mined by u>x , and we obtain a dimensionless parameter 
characterizing the relaxation: 

(T cWl MVa)(VAM 3 / 5 (oV^i) 2 . (23) 

This estimate shows the possibility of both fast and slow 
relaxation. For t c u>i 3> 1 one can expect to observe osci- 
ullations of the conensate. 

For a small interaction (77 ps Nq<x/Iq <C 1) the charac- 
teristic energy of expansion ~ Hloq , and we have {t c ui\ )~ 
(Iq / o){Iq / Nqo)(u)q / LO\) 2 3> 1. In this limit the relaxation, 
and hence the loss of coherence, is always slower. 
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